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1. Introduction 

The question of how to gain information about the state of a quantum system has intrigued 
researchers since the early days of quantum mechanics. With the advent of quantum 
computation, this fundamental question became also of practical interest, mainly because 
the final stage of a quantum algorithm necessarily is qubit readout. This task only requires 
distinguishing between two particular qubit states and, thus, can be achieved by projective 
measurements. Nevertheless, going beyond readout is of interest as well, since one 
also desires direct experimental evidence for coherent superpositions emerging, e.g., from 
tunnelling oscillations. 

In order to obtain a quantum mechanical description of a measurement process, one 
usually models the measurement apparatus as a macroscopic quantum environment, i. e., as 
a heat bath, where the pointer of the apparatus corresponds to an effective bath coordinate. 
When interacting with the central quantum system, the bath acquires information about the 
system state. Owing to the macroscopic nature of the bath, one may assume that already a 
fraction of the bath possesses the full information about the effective pointer coordinate yj. 
Therefore one can obtain knowledge of the pointer position without violating fundamental 
laws of quantum mechanics. 

Recently, superconducting quantum circuits have provided a new arena to test 
fundamental questions of quantum mechanics in the laboratory. Prominent examples are the 
demonstration of coherent time evolution in charge qubits |2J and of Berry phases Q, as well 
as testing Bell inequalities ID. Above all, different protocols for quantum measurement were 
successfully implemented in circuit quantum electrodynamics ||5}{8l. For a superconducting 
solid-state qubit, the practical measurement of one of its coordinates is performed by 
coupling it to a macroscopic environment, as given by external circuitry, via a quantum point 
contact Il9l[l0l or a harmonic oscillator Depending on the setup, the oscillator is realised 
by a dc superconducting quantum interference device (SQUID) ifTTl or a superconducting 
resonator ifTZl . In both cases, the resonance frequency of the oscillator depends on the 
qubit state. Consequently, the response of the oscillator to a close-to-resonant ac-excitation 
possesses a phase shift which can be measured, and from which one can infer the qubit state. 
First experiments in this direction worked with an oscillator whose frequency was much lower 
than the qubit splitting ||5][l3l[T4l . More recent experiments ||6l[T5l operated in the so-called 
dispersive regime, where the oscillator frequency and the qubit splitting are of the same order, 
while their detuning is still larger than their mutual coupling strength. A crucial detail is 
that the oscillator frequency and bandwidth naturally limit the time resolution in such qubit 
measurements. Thus, using the said schemes with slow oscillators, it is only possible to 
extract time-averaged information about the qubit state in general, but there is no possibility 
to resolve its dynamics in time. 

Recently, a first step towards a time-resolved measurement of qubit dynamics has been 
proposed II16II17I : When a weak high-frequency field acts directly on the qubit, the reflected 
signal acquires a time-dependent phase shift by harmonic mixing. Lock-in amplification of the 
reflected signal then allows obtaining information about the qubit dynamics. In this work, we 
combine both approaches and extend the scheme of Refs. 01611711 to a qubit coupled to a driven 
high-frequency oscillator A measurement protocol for such a setup is particularly appealing 
because an oscillation mode is part of most recent superconducting qubit designs. Moreover, 
the oscillator serves as filter for quantum noise and, thus, reduces qubit decoherence. Here, we 
focus on a flux qubit embraced by a dc-SQUID, whose fundamental frequency may even be 
tunable to some extent 11811191 . As a particular feature of this realisation, the qubit-oscillator 
coupling is non-linear in the oscillator coordinate, that is, the coupling possesses both a 
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Figure 1. Sketch of the flux qubit (blue) coupled to a dc-SQUID. The interaction is characterised by 
the linear coupling g\, which depends linearly on the SQUID bias current /b, and the quadratic coupling 
g2- The SQUID with Josephson inductance Lj is shunted by a capacitance C. The frequency shift of the 
resulting harmonic oscillator (green) can be probed by external resonant ac-excitation Acos(tla(;f) via the 
transmission line (black), in which the quantum fluctuations £,^^{t) are also present. 



significant linear and quadratic contribution. It will turn out that for realistic parameters, 
the measurement scheme relies on the quadratic part of the coupling. 

The paper is structured as follows. In section |2] we introduce our model and 
discuss dispersive qubit readout in generalised terms. The central relation upon which our 
measurement scheme relies is derived in section|3] Section|4]is devoted to numerical studies 
in which we test our measurement relation and work out quantitatively measurement fidelity 
and backaction. Furthermore, we provide an estimate of the signal-to-noise ratio and discuss 
realistic parameters for a possible experimental implementation. The appendix contains 
details about the derivation of the measurement relation, the input-output formalism ll20l and 
the Bloch-Redfield master equation which we use for obtaining numerical results. 



2. Dissipative qubit-oscillator model 

2.1. System-bath model 

We consider a superconducting flux qubit coupled to a SQUID |l7| as sketched in figure[T] The 
SQUID is modelled as a harmonic oscillator, which gives rise to the Hamiltonian l6i.7l [T8ll2TI 

^ ^ ^^(j^^^jj^l^^t^^ 1^ +r7(c7-cos0-o:,sin0) [gi{a + a'^)+g2{a + a'^f] . (1) 

The first term represents the qubit with energy splitting TiOqb = fi{e^ + 5^)'/^ and the mixing 
angle 9 = arctan(5/e) which depends on the controllable qubit bias energy e and the qubit 
gap energy 5, while Ox.z denote the Pauli matrices. The second term describes the oscillator 
with frequency Q. and the bosonic creation and annihilation operators a^ and a, respectively. 
The qubit couples to the oscillator in two ways. First, via dipole interaction with strength gi, 
which is linear in the oscillator coordinate a + fl^. Up to order g^, this causes a frequency 
shift for both the oscillator and the qubit. The second coupling term proportional to g2, by 
contrast, is quadratic in the oscillator coordinate. Its physical origin is a non-linear Josephson 
inductance which depends on the magnetic flux, by which the SQUID is penetrated iTTl fTSl . 
This term provides a frequency shift already in first order of g2. The interaction coefficients 
gi and g2 can be controlled to some extent, as an expansion of the qubit-SQUID interaction 
to second order in the oscillator coordinate demonstrates fflSl . In detail, for a small SQUID 
bias current 4, the coupling coefficient gi depends linearly on /b, whereas g2 is independent 
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of the latter, as we discuss in section |4!4l For 4 = 0, one even obtains g\~Q, such that the 
qubit couples only to the square of the oscillator position llTl lTSl . 

Regarding a time-resolved measurement of the qubit dynamics via the oscillator, it will 
turn out that for realistic parameters of flux qubits, this quadratic coupling is crucial, while 
the linear coupling turns out to be typically too weak. For common circuit-QED setups using 
charge and flux qubits coupled to a transmission line resonator 11 Illl2ll22ll23l , not only but 
also g2 is too small. Thus, we henceforth focus on setups of flux qubits coupled to SQUIDs 
possessing a sizable quadratic coupling, as described above. 

The qubit-SQUlD system is further coupled to external circuitry, which acts as a 
dissipative environment and is modelled by the system-bath Hamiltonian Il24l - l27l 




Here, Q = a* + a is the oscillator coordinate, such that the interaction term represents 
the inductive coupling between the qubit and the flux degree of freedom of the SQUID. 
The system-bath interaction can be fully characterised by the spectral density J{(o) = 
KY,k kA:p5(a) — Ojt) which is proportional to the real part of the effective impedance of the 
environment 1281 . Here we assume an ohmic spectral density, 7(ft)) = ao), for which the 
dimensionless damping strength a can be interpreted as effective resistance ||29l - l3n . 



2.2. Qubit-oscillator interaction in the dispersive limit 

We are interested in the dispersive limit which is characterised by a detuning A = il — tOqb 
larger than the qubit-oscillator couplings, 

gi,g2<|A|, A = i2-«qb. (3) 
It is then convenient to go to the dispersive picture via the unitary transformation (IA.2) . As 



detailed in Appendix A this yields the effective Hamiltonian 

Jfo = '^''Jfo'^=h£l(d'd+^ + ^^^o-, (4) 

where the transformed bosonic operators a and are defined in equation dA.lOb . The qubit- 
oscillator coupling has been removed formally by shifting it to the operator-valued oscillator 
frequency 



40) 

n = i2yi + — , (5) 
where the overbar denotes the dispersive picture, while the qubit operator 

(b =^a-(- i ^ sin^e + ^aJ- i ) cos sin 



2 -VA il + OJqb/ 2 VA i2 + fi)qb/ (6) 

+ g2{(y-cos,Q - (7jSin0), 

determines the coupling. The interpretation of equations (|5]l and ^ is that the oscillator 
frequency depends on the qubit state. This allows dispersive qubit readout by measuring the 
associated phase shift of the oscillator response upon resonant driving, like in the case of a 
classical oscillator that is driven by an external force. In particular, assuming cos 0=0 and 
g2 =0, equation (|5l) predicts the frequency shift A = £2 + a;§j [1/A- 1 / (ii + ojqb)]^'. The last 
contribution in £2 stems from counter-rotating terms in the qubit-oscillator interaction. These 
must be accounted for the case of large detuning A where a rotating-wave approximation 
produces inaccurate results ll34l . Depending on the qubit expectation value (c-), the oscillator 
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is red or blue detuned. Thus, we obtain in this Hmit the well-known qubit-dependent phase 
shift corroborated in various experimental realisations [I5l 48l22l . There, however, the oscillator 
frequency was smaller than the qubit splitting, i2 ^ Oqb- As a consequence, it was only 
possible to obtain time-averaged information about the qubit state. 

Now the goal of this paper is a generalisation of dispersive qubit readout such that 
time-resolved information about the qubit state can be obtained as well. This obviously 
requires oscillator frequencies and bandwidths larger than the qubit transition frequency, that 
is, Q. ^ (Bqb. We emphasise that equations (|5]) and (|6]l are nevertheless valid as long as the 
coupling constants are small enough to fulfil condition (O; for details see Appendix A and 
Ref. [34|. If the qubit dynamics is much slower than the oscillator, the qubit can be treated 
within an adiabatic approximation. This means that the qubit dynamics is assumed to be 
constant during one oscillator period. In turn, the time evolution of the oscillator depends on 
the instantaneous qubit state. Then the Schrodinger-picture operators Ux^z in equation (|6]l can 
be replaced by their time-dependent expectation values, and the operator- valued quantity d) is 
substituted by 

, , (§isin0)2 / I 1 \ , , /I 1 A , ^ 
(b{t) = ' - (c7,), + 2i cos sine (a,) 



2 \A Q. + COqbJ 2 \A a + COqb/ (7) 

+ g2(cos0(c7,), -sin0(c7,),) . 

Equation (|7]) implies that information about the time-dependent qubit state is encoded in the 
effective oscillator frequency Q. = n(f ). This gets expressed as a slow parametric modulation 
in time, Uke for a parametric oscillator In detail, the instantaneous qubit state enters via the 
qubit expectation values ((J,,,-)? = Tiq\,{ax,zPQ{t)}, where Trqb denotes the partial trace over 
the qubit degrees of freedom. The time dependence, indicated by the subscript (. . .),, stems 
from the evolution of the total qubit-oscillator state po(f) under the effective system-bath 
Hamiltonian ( IA.15I ). 

As an important intermediate result, the found modulation of Q. in time can be traced 
back to the qubit dynamics in the absence of the driving. This enables measuring the qubit's 
time evolution via the oscillator response to resonant driving. 



3. Time-resolved measurement of the qubit dynamics 

The qubit-oscillator Hamiltonian in the dispersive picture, equation (|4|, together with the 
effective, modulated frequency (|5]l already indicates that the oscillator detuning may contain 
information about the qubit dynamics. As in the case of the traditional dispersive readout, 
we consider the response of the system to an ac-field that is resonant with the oscillator 
Physically, the situation corresponds to that of a classical mechanical oscillator driven by 
an external periodic force. Owing to the only weak dissipation, the response is manifest in 
the phase of the reflected ac driving. In the following, we establish a relation between this 
phase and a time-dependent qubit expectation value. This relation will form the basis of our 
measurement protocol. 



3.1. Response of the qubit-oscillator compound to resonant driving 

In the theory of optical cavities, the response to an external ac excitation is conveniently 
calculated with the input-output formalism 12011351 . This formalism has also been appUed to 
quantum circuits ||5][l3][l4l|T6l . Its cornerstone is the relation 

^out(f)-^in(0 = 2aG, (8) 
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formulated in the Heisenberg picture and derived in [Appendix B[ It relates the incoming 
and the outgoing fluctuations of the transmission line, ^in/out(0' the time-derivative of the 
system-bath coupling operator, which in our case \s Q~a + a\ The dimensionless dissipation 
strength a of the ohmic spectral density quantifies the coupling between the oscillator and 
the electric environment and, thus, appears as prefactor An ac -driving corresponds to a 
coherently excited incoming mode, such that the fluctuations can be separated into quantum 
fluctuations ^^™{t) and a deterministic component A cos (iiacf)- Here, the deterministic part 
is an ac-field in resonance with the bare oscillator, i2ac = ii, such that 

^in(0 = C(0+Acos(nf), (9) 

which implies the expectation value (^in(f)) = Acos(nf ). Then the input-output relation (O 
becomes <^oiit(0 = ^^™{t) +Acos(i2f) +2aQ and, thus, the corresponding expectation value 
of the outgoing signal reads 

(<^out(0)-^cos(i2f)+2a(e). (10) 

Also here, it is convenient to work in the dispersive picture obtained by the unitary 
transformation (IA.2t . While this leaves the environment operators unchanged, the coordinate 
by which the oscillator couples to the environment changes as Q^Q = d + a'— (Aa — 
'^i)<7i +2AnCr.; see equation ( IA.14b . The time-derivative dg/df can be obtained from the 
commutator of Q with the Hamiltonian (01) augmented by a term that describes the driving. 
This yields terms of the order i2 and terms with prefactors OJqb and g\/A. For a fast oscillator, 
the latter terms can be neglected, and we obtain the equation of motion 

^e + 2an^e + n2e==-2n[^i^„'"(f)+Acos(nr)] . (11) 

This linear, inhomogeneous equation is readily solved with the help of the Green's function 
for the dissipative harmonic oscillator Inserting the resulting dQ/dt into the input-output 
relation ( fTOl i and neglecting transient terms yields for the expectation value of the outgoing 
signal the expression 

(^out(0)=^cos{i2f-(p(f)}. (12) 

Owing to the weak dissipation, the system energy is almost preserved, such that the amplitude 
of the incoming and the outgoing signal are practically the same. The phase difference 

stems from the coupling to the qubit which detunes the oscillator, while the slow time 
evolution of the qubit renders the phase time-dependent. The approximation is valid if the 
qubit-oscillator couplings are smaller than the oscillator damping rate, i.e., for gi,g2 <C CcQ.. 
In other words, the first term in the denominator is negligible, since the qubit-induced 
frequency shift H — H is of the order gi2- This also ensures (p I and, thus, (p w tan 9. 
Next, we insert the effective frequency (|5) together with equation d?) and obtain to second 
order in gi and first order in g2 the phase shift 



2^? ^ 1 1 



"^^'^ - ^ [l-^^J [-n^e(a,),+cos0sin0(a.), 
+ ^(cos0((J,),-sin0((J,),) . 



(14) 



This central relation forms the basis for our non-invasive qubit measurement via a resonantly 
driven harmonic oscillator It identifies a set of qubit observables, which generate the low- 
frequency system dynamics, as the cause of a small phase shift between the ingoing and 
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outgoing signal. In other words, equation ([14} enables one to monitor the qubit dynamics by 
continuously measuring the phase shift (p(t) with suitable experimental techniques. 

By evaluating the prefactor for specific setups, we will see below that our measurement 
scheme is particularly feasible for flux qubits. In this case, the last term of the phase shift ( fT4] i 
dominates, and one measures the qubit variable <T^cos0 — CTvSinS, i.e., the flux degree of 
freedom by which the qubit couples to the SQUID; cf. the model Hamiltonian ([T|). 



3.2. Static versus dynamical phase shift 

The terms entering the phase shift <p{t) may be static as well as dynamical. In the first instance, 
this depends on whether or not the related qubit observables undergo any time evolution. At 
this point, further insight is obtained by a closer look to the Heisenberg equations of motion 
for the qubit operators (Jx and cr^. In the dispersive picture, valid under conditions (|3}, they 
are derived from the effective Hamiltonian given by equation (IH, and read 

Ox = ^[^,0:J = -ajqbCTy , (15) 

(T, = ^[^^o,ctJ=0. (16) 
h 



Thus, in the dispersive qubit-oscillator coupling limit (see section 12. 2t . the observable C7, 
is a constant of motion. As a consequence, those contributions to (p{t) that depend on 
((7;)r = (o^;)const are time-independent. This corresponds again to the established scheme 
for non-invasive qubit state readout. 

On the contrary, the observable Ox possesses a non-trivial time dependence generated by 
Thus, {ax)t renders the phase shift (p{t) dynamical. This, in turn, enables a time-resolved 
single-run measurement of the unitary qubit evolution by means of the qubit observable Ox- 
According to our measurement relation (Il4l . the dynamical phase signal has the amplitude 



2 

•Pma 



an 



??sin0cos0 f I — 2?2sin0 



(17) 



Interestingly, <p^^^ is reciprocal to the damping strength and the oscillator half-bandwidth 
aii. Thus, a large oscillator frequency together with strong damping lead to reduced angular 
visibility. On the other hand, the adiabatic treatment of the qubit underlying relation (fl4] i 
becomes invalid if either i2 or a are too small. Moreover, the input-output relation (ITO} 
crucially relies on finite damping. Thus, appropriate choices for £1 and a need to be based 
upon a compromise between good phase resolution and the validity of our approximations. 
We go into further detail about this issue when discussing the measurement quality in 
section |4!2] 

It is important to note that the amplitude ^^^^ possesses contributions from both the 
linear and the quadratic qubit-oscillator interaction of Hamiltonian ([T]). The first term on 
the r h. s. of equation ( fTTI i stems from the linear interaction characterised by the coupling 
coefficient g\. Like the effective Hamiltonian ( IA.13I ), this contribution is of second order 
in the dispersive parameter gi/A. Due to the minus sign inside the round brackets, it is 
additionally minimised, given that A/(I2 + Oqb) « 1 for large detuning A. Thus, for a Cooper- 
pair box or a flux qubit coupled to a high-frequency transmission line resonator, where the 
qubit-oscillator coupling merely depends on the circuit characteristics and therefore is purely 
linear, i. e. ^2 = 0, the maximum amplitude (p^^^ drops below any useful level. 

On the contrary, a finite quadratic qubit-oscillator interaction g2>0 ensures a noticeable 
phase signal, independent of the detuning A, and even if g2 ^ ^i- Formally, this feature 
arises from Hamiltonian where the quadratic coupling term already generates a frequency 
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shift in zeroth order perturbation theory. If the qubit-osc illator interaction is transverse, that is, 
£ = 0, the phase resolution is maximised, whereas it vanishes for purely longitudinal coupling, 
i.e., for 5 = 0. Hence the presence of a non-linear qubit-oscillator interaction, as provided by 
a non-linear SQUID Josephson inductance, turns out to be a crucial ingredient. 



4. Measurement quality 

Still it remains to corroborate the central measurement relations ( fT4l i and ( [TtT i. respectively, 
via comparing them with the phase shift obtained by simulating the actual measurement 
process. In doing so, we restrict ourselves to the fundamental example of coherent qubit 
oscillations. For the numerical treatment of the qubit-oscillator state, we employ the quantum 
master equation ( IC.lb derived from the full dissipative qubit-oscillator-bath Hamiltonian 
For a realistic evaluation, we use parameters similar to those of the experiment reported in 
references II18II191 . Furthermore, we employ A' = 10 oscillator states, which turns out to be 
sufficient to reach numerical convergence. 



4.1. Time-resolved measurement of unitary qubit evolution 

If the qubit is only weakly coupled to the oscillator, and if the latter is driven only weakly, 
the qubit's time evolution is rather coherent (see section |43] on qubit decoherence). For this 
scenario, figure |2j a) depicts the time-dependent phase ^(f) computed with the measurement 
relation (fl4] i. while the inset confirms its proportionality to the qubit expectation value {Ox)t- 
For a comparison, we wish to recover this phase information directly by analysing the 
outgoing signal (^out(O)' given by equation ( fTOl i. In an experiment, this can be achieved 
by lock-in techniques which we mimic in the following way lf36l : First, we focus on the 
associated spectrum ((^out(fi^)) depicted in figure|2jb). It reflects the qubit dynamics in terms 
of two sidebands around the central peak related to the oscillator frequency, here chosen 
as £2 = lOOqb. The dissipative influence of the environment, modelled by a transmission 
line (see figure [U, is reflected in a broadening of this peak. The corresponding oscillator 
bandwidth is given as 2a£2, where a denotes the dimensionless damping strength; see 



Appendix C Here, we recall that the oscillator is driven resonantly by the external driving 
signal Acos(£2acf)' that is, £2 = £2ac. In the time domain, the sidebands correspond to the 
phase-shifted signal (^out(O) = Acos(£2f — (Pexp(O) with slowly time-dependent phase 9exp(f )■ 
In order to obtain this phase (Pexp(0' we select the spectral data from a frequency window of 
size 2A£2 centred at the oscillator frequency £2, which means that ((^out(o)) is multiplied 
with a Gaussian window function exp(— (o — £2)^/A£2^). We choose for the window size the 
resonator bandwidth, A£2 = a£2, which turns out to suppress disturbing contributions from the 
low-frequency qubit dynamics. Finally, we centre the clipped spectrum at zero frequency and 
perform an inverse Fourier transform to the time domain. If the phase shift 0ext was constant, 
one could use a much smaller measurement bandwidth. Then the outcome of the measurement 
procedure would correspond to homodyne detection ll35l of a quadrature defined by the phase 
shift and yield a value oc cos ^exp- 

Figure |2j a) reveals the good agreement of the resulting (pexp(0 with the prediction of 
our measurement relation, (p(f ) oc (cTv),, at angular resolutions of 1-2°. Good agreement is 
already obtained for a oscillator frequency £2 = lOcOqb, which obviously represents a good 
compromise between the validity of the adiabatic approximation (see section I2.2l i and a 
sufficiently strong signal. There is even some room for obtaining a stronger phase signal 
since the dissipation strength a still can be reduced without violating the validity range of our 
theory as long as Oqb < A£2. 
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Figure 2. Time-resolved measurement of coherent qubit oscillations at the degeneracy point e = 0. The 
full qubit-oscillator state was simulated with the quantum master equation l IC.U with A' = 10 oscillator 
states and the parameters £1 = il^c = lOtOqb, gi =0.1 Oqti, g2 = 0.01 Oqti, A = l.OcOqb. The dimensionless 
oscillator dissipation strength is o; = 0.12. The resonator bandwidth is given by 2aQ. = 2.4£0qb. (a) Lock- 
in amplified phase <Pcxp{f) (dashed green hnes), compared to the estimated phase (p(t) (solid red line) of 
the outgoing signal {i^out('))- Here, (fi(t) « (cTv)r [cf. equation )14H . which is corroborated by the inset 
showing that (fTv)/ performs oscillations with (angular) frequency co^^y. (b) Power spectrum (i5oui('h)) for the 
resonantly diiven oscillator (blue solid line). The sidebands stemming from the qubit dynamics are visible 
at frequencies £2 ± (Bqb. In order to extract the phase information, we apply a Gaussian window function 
with respect to the frequency window of half-width Ail = 1.2cBqb, which turns out to be the optimal value 
for the measurement bandwidth. 



Setting either g\ or g2 to zero (not shown) reveals that the non-linear coupling g2 is 
responsible for the good agreement of the phase shifts in figure|2ja). Thus, the whole protocol 
is mainly applicable to flux qubits coupled to SQUlDs. For charge qubits, by contrast, the 
typical values for g2 are too small. Furthermore, we have verified that the visible constant 
delay between both phases (p[t) and 9exp(0 does not depend on the selected parameters, 
while its detailed origin remains unexplained. 

4.2. Measurement characterisation: fidelity and backaction 

The validity of relation (fl4l i for the phase 9(f) is naturally limited to specific parameter ranges 
due to the various underlying approximations made. The main crucial assumptions to justify 
the adiabatic treatment of the qubit are a large qubit-oscillator detuning A/O ~ 1 and weak 
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Figure 3. (a) Fidelity defect 8F= 1 —F for the phases (f>(t) and <Pcxp(0 and (b) time-averaged trace distance 
D between the density operators of a qubit with finite coupling to the oscillator and a reference qubit without 
oscillator. Both quantities are depicted for various coupling strengths gi in dependence of the oscillator 
frequency Q.. All other parameters are as in figure |2] 



mutual interaction, gi A. Furthermore, the oscillator damping a is assumed to stay 
within the limits gi2/n<ca-Cl. 

In an experiment, the oscillator frequency and the coupling strength are finite, though. 
Consequently, the actual phase ^exp(0' which we extract numerically and which can be 
measured by lock-in amplification, generally differs from the predicted phase <p{t). Thus, 
the mutual agreement of both phases needs to be tested quantitatively for realistic scenarios. 
To this end, we employ the measurement fidelity F with the scaled overlap defined as 

1-1/2 



(?',9'exp) 



Jdt(p'^(t) Jdt(p^ 



(0 



d?<?)(0<?>exp(0 



(18) 



The ideal value of F = 1 is assumed if (p{t) oc (p^^^(t). 

In figure [3 a) we depict the fidelity defect 5F = l—F between (Pexp(0 and ^(f) as 
a function of the oscillator frequency Q. = iiac for different quadratic coupling coefficients 
g2- As expected, the overall fidelity is rather insufficient for small oscillator frequency 
Q. < lOOqb, for which the adiabatic approximation of section l272] is not valid and, moreover, if 
the oscillator bandwidth is too small to resolve the qubit dynamics, i.e., if OJqb < 0Ci2. Along 
with increasing Q., the fidelity defect 8F drops to values of 0.05-0.5, independently of the 
parameter g2- Taking into account that the fidelity is arbitrarily lowered by the constant delay 
between (p{t) and (Pexp{t) visible in figure|2] this still corroborates that Q. — lOcOqb is a good 
choice. In the Umit of large oscillator frequencies, we again observe an increase of the fidelity 
defect, which occurs the sooner the smaller g2- This latter effect, which is only visible for 
the smallest value of g2 in figure [3ja), is directly explained by a reduced maximum angular 
visibility of the phase (p{t) oc gj/Q.. Thus, figure |3j a) provides a pertinent indication for the 
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validity frame of our central relation (IT4l . 

Moreover, it is necessary to take into account the backaction upon the qubit that stems 
from the non-linear qubit-oscillator interaction. An appropriate measure for how much the 
qubit dynamics is perturbed by the oscillator is given by the time-average D of the trace 
distance D{t) = 7Tr|pqb(f) — Pqb,o(OI between the qubit dynamics with and without the 
coupling to the driven oscillator To be specific, we compare the qubit state Pqb(0 evolving 
under the full system-bath Hamiltonian (|2]i to an unperturbed reference state Pqb,o(0 that 
evolves unitarily under the bare qubit Hamiltonian ,i^b = (^2t'^qb/2)c7-. Thus, the trace 
distance essentially quantifies the invasiveness of the measurement based upon the second- 
order qubit-oscillator interaction. In the absence of perturbations to the qubit, D vanishes by 
definition, while D = 1 if the density operator of the measured qubit is completely unrelated 
to that of the reference. 

Figure[3jb) shows that the predicted phase ^(f) faithfully describes the unperturbed qubit 
dynamics as long as the coefficient g2 stays sufficiently small. A reliable operating range 
appears to be g2/(i>qh ^ 0.01. For £1 = 10 ©qb, this is consistent with our above reasoning 
regarding the fidelity. For even weaker quadratic interactions, we first find D oc which 
implies that the dispersive first-order coupling in terms of gi governs the qubit-oscillator 
interaction when Q. is small. This cubic dependence is due to relation ( fT4l i and to the fact 
that the first-order perturbation acting on the qubit has an inverse quadratic dependence on 
the detuning A oc f2. Beyond a critical detuning, which individually depends on gi, the 
quadratic interaction prevails again, as is reflected by the saturation of D with increasing 
oscillator frequency £2. For g2/(0qb ^ 0.05, the effect of the linear coupling withers. Thus, at 
g2 = 0.01 ©qh, the value of gi is rather irrelevant for our measurement protocol. 



4.3. Signal-to-noise ratio 

Generally, a signal can be resolved only if its spectral density exceeds the level of background 
noise at the measurement frequency. In the present scheme, the desired information is 
contained in the sidepeaks of the spectrum at £2 ± ©qb; see figure|2tb). For the corresponding 
phase-modulated oscillation Acos{Q.t — (Pmax sin(©qbf))i these sidepeaks correspond in the 
time domain to the oscillation A(j!)niaxsin(©qbf). If the signal is integrated for a time t, its 
spectral weight becomes (A^niax)^f - The phase amplitude ^max is given by equation (T71 . but 
for the present purpose, it is sufficient to consider the dominating contribution which is the 
one proportional to g2- Thus, here ^niax = while we restrict ourselves to the case 

e = 71 /2. 

Since the measured signal corresponds to the state of a highly excited environmental 
mode, the relevant noise level is determined by the fluctuations of the effective bath coordinate 
^ . If the temperature is sufficiently low, such that thermal excitations do not play any role, its 
spectral density equals the bath spectral density: S^^{co) = /(©) = a©. Thus, the signal is 
at least as big as the noise background if /(£2) < (A(jciinax)^f- In other words, the time t during 
which the output is recorded must fulfil 

47(n) _ (an)3 _ 

-(A<p,.ax)2~(2^2A)2-'™^^- ^'^^ 

For the parameters used in figure |2] f^eas ~ 4 • 10^/ ©qb. 

Since the measurement is performed via a coupling to external degrees of freedom, the 
qubit experiences unavoidable decoherence, which means that coherent qubit oscillations fade 
away with a decoherence rate 7^. The time during which meaningful information can be 
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obtained is therefore limited by the inequality f < 1 /y,/, . This condition together with condition 
(fT9l ) can be fulfilled only if y^Was > 1- 

The decoherence rate can be estimated upon noticing that our qubit Hamiltonian 
represents a generalised spin-boson model [24-26] with the bath coupling ^cTiTJ, where 
Tj = IgiQ^ is the effective bath coordinate. For weak dissipation, is given by the auto- 
correlation function of the latter evaluated at the qubit splitting, i.e., jif, = Crm (Cqb) (371. This 
still holds true in the presence of ac driving provided that the driving-induced renormalization 
of the qubit splitting is negligible f38|. We separate the qubit coordinate into the responses 
to the driving and to the incoming fluctuations, Q ~ (Q) + 5Q; cf. equation Q. Then we 
proceed along the lines of Ref. Ii39j : The relevant terms are those of second order in 5Q, such 
that Crjriit) ~ 4g2{A/aD.)'^cos{Q.t){5Q{t)5Q{0)). By Fourier transformation, we obtain to 
lowest order in (Oqb / ^ the decoherence rate 

n (20) 

which is the inverse of the required measurement time fmeas- A comparison with the 
numerically computed decay of the qubit coherence (not shown) confirms this value. The 
obtained relation y^fmeas = 1 marks the quantum limit of a measurement [391 and allows one 
to marginally fulfil both conditions on the measurement time t. 

4.4. A possible experimental implementation 

Specific parameters can be obtained for the setup of reference |fT9l for which the qubit- 
oscillator coupling parameters are determined by the flux bias current I\, and read 

M/p sin((p/2) 

^' 4S/ccos2((p/2)Y2Lj ^ ' 

MIp sin(^/2) 
16Lj/ccos2((p/2)" 

Here, M denotes the mutual SQUID-qubit inductance, /p is the qubit persistent-current, Iq 
the critical current of the SQUID Josephsonjunctions, and Lj = 0o[47r/ccos(;r0sQ/0o)] ' the 
SQUID Josephson inductance. The flux ^sQ that penetrates the SQUID loop corresponds to 
the phase <p = 2k^^q/^q, with 0o = /i/2e being the flux quantum. For small bias currents 
4 while neglecting the inductance of the wire that leads to the shunting capacitance C, the 
oscillator frequency is approximately given by the SQUID plasma frequency Q. = jLjCj^'/^ 

For a typical qubit transition frequency of tOqb /2k = 5GHz, resolving the qubit dynamics 
requires an oscillator frequency of £l/2n ~ 10a)qb/2;r = 50GHz. The necessary lock-in 
measurements at a carrier frequency £l/2n ~ 50GHz are particularly challenging at low 
temperatures, but feasible BOl . They require rather expensive amplifying technology such as 
cryogenic amplifiers developed by Low Noise Factory (Sweden). Recently developed suitable 
devices can be operated at up to 36GHz and possess reasonably low noise temperatures. 

Alternatively, a Josephson parametric amplifier 11411421 enables the detection of oscillator 
frequencies as high as Q./2K ~ 20-25GHz. Thus, as a compromise, we restrict ourselves 
to an oscillator frequency of Q./2n = 24GHz, which is suitable to detect the dynamics of 
a qubit with a)qb/27r = 2.5GHz, a value still large enough to avoid thermal excitations at 
working temperatures of 20mK. The above value for £1 can be realised using the parameters 
/c = 4.25/iA, <p = 2.37r, and C = IpF, which are similar to those of references [1T8][T91 . 
Along with M ~ 17.5pH, /p = 300nA and 4 = 0.4/xA, we obtain the qubit-oscillator coupling 
coefficients g\l2n — 45MHz and g2 — 24MHz. Thus, the relevant dimensionless couphng 



^2 = -TT7^— I^ii- (22) 
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assumes the value ^2/ Cqb — 0.01 used in our numerical studies. An adequate oscillator half- 
bandwidth is An = 2.9GHz, which implies a low external quality factor of 2 — 4-5. 

Our two-state model for the qubit does not consider possible excitations to non-qubit 
states caused by the coupling to the oscillator Nevertheless our modelling is appropriate, 
because such leakage has far less relevance than for a Cooper-pair box, owing to the fact 
that the higher states couple only weakly to the SQUID B3l . Apart from this, it is possible 
to design or tune the oscillator such that its frequency is far from any qubit resonance. The 
required oscillator frequency of the order lOGHz is still significantly smaller than the gap 
energy of aluminium, such that quasi particle excitation should not play a major role. This 
issue is even less critical for niobium. 

5. Conclusions 

We have generalised dispersive qubit readout to time-resolved observation of the qubit 
dynamics. Concerning the setup, the main difference to dispersive readout is that in the 
present proposal, the oscillator frequency needs to exceed the qubit splitting by roughly one 
order of magnitude, and the oscillator bandwidth should be at least twice the qubit frequency. 
Also here, the oscillator frequency becomes dynamically red or blue detuned, depending on 
the state of the qubit. When driving the SQUID oscillator at its bare frequency i2, this detuning 
turns into a phase shift visible in the reflected signal via lock-in techniques. For such qubit 
measurement using the oscillator phase, the oscillator frequency represents the sampling rate, 
which explains the need for high frequencies. 

The constituting measurement relation has been derived from the input-output formalism 
under time-scale separation of the bare qubit dynamics from the oscillator. A numerical 
solution of the Bloch-Redfield master equation for the full qubit-oscillator dynamics allowed 
us to compute the phase of reflected signal also directly. Its good agreement with the phase 
predicted by our measurement relation confirms the validity of the latter even when the 
oscillator frequency is just moderately large. Thus, there is no need for driving the qubit 
with extremely high frequencies, which would be quite challenging in an experiment. The 
found agreement is also reflected by the measurement fidelity, which already for moderate 
frequencies is rather good. Furthermore, the numerical analysis has demonstrated that the 
external ac-driving does not significantly modify the qubit dynamics, which means that the 
backaction of the measurement process is weak. However, it must be emphasised that the 
whole scheme relies on the coupling of the qubit via the oscillator to a dissipative environment, 
which causes qubit decoherence already when the external driving is not active. In the limit 
of far qubit-oscillator detuning, this qubit decoherence gets drastically reduced though. 

Evaluating the measurement relation for parameters of recent experiments with flux 
qubits predicts phase shifts up to 2°, which can be measured. Moreover, it reveals that the 
signal mainly stems from the coupling of the qubit to the square of the oscillator coordinate. 
The linear coupling to the coordinate, by contrast, leads to a rather small phase shift. This 
means that the measured quantity is essentially the qubit's flux degree of freedom. Likewise, 
the linear coupling of a superconducting charge qubit to a waveguide resonator is also too 
weak. Since for this system the non-linear coupling practically vanishes, the measured signal 
remains tiny. In conclusion, with present technologies, our measurement protocol should 
be feasible best with flux qubits coupled to SQUlDs that possess a significant non-linear 
Josephson inductance. All in all, our proposal may initiate further progress on the way towards 
single-shot experiments that demonstrate quantum coherence in solid-state devices. 
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Appendix A. System-bath Hamiltonian in the dispersive coupUng Umit 

In the limit of large oscillator-qubit detuning, the coupling coefficients automatically fulfil the 
conditions 

gl,g2<A, A^il-OJqb, (A.l) 

which mark the dispersive coupling regime. Following references ll32ll33l . the effective 
Hamiltonian ,i^,disp = is then obtained from the full system Hamiltonian ([T]i by 

the unitary transformation 

^ = exp (Aa^ + AiJ5^ + XaW) , (A.2) 

where 

^ = a^fl''' - c7+fl , (A.3) 

,y = a- a - , (A.4) 

r = (7-(fl-fl^). (A.5) 

Defining Z = Oqb + i2, the necessarily small and dimensionless dispersive parameters 
gisine 

Aa = — , (A.6) 

A 

?isin0 

Ae= (A.7) 
,?i cos Q 

Aa = -^ (A.8) 

emerge. Expanding the transformed Hamiltonian in powers of AAx.n. we obtain to second 
dispersive order the effective Hamiltonian 



^ = hO. [a^a + 2 j + + 2 (^^A - ^A|) 0-{a + a 

+ ^Aa (Aa + Ar) (Tv(« + a'f- ^Aa (Aa + Ar) ay(a' - {a'^f) (^-9) 

+ g2 (cos CT, — sin C7v) (fl + fl^) ^ • 

The third and fourth term of this Hamiltonian constitute corrections to the curvature of the 
oscillator potential, i.e., the prefactor of [a + a^)^. They stem from the linear qubit-oscillator 
interaction and, thus, enter only in second dispersive order Since we consider a high- 
frequency oscillator, the detuning is always positive, A > 0, such that the dispersive parameters 
Aa and Aj; are of opposite sign. Thus, in the case of far dispersive detuning £l ^ a)qb, these 
terms become rather small. In spite of this, we keep them up for later purpose. On the contrary, 
we can safely neglect the fifth term which is not of the shape {a + a^)^ and whose coefficient 
is small as compared to the other terms. 
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The last term of equation (|A.9t . stemming from the second-order interaction between 
the qubit and the oscillator, plays a particular role. Since it already is of second order in 
the oscillator coordinate a + and its coefficient g2 is correspondingly small, g2 g\, it is 
not affected by the transform (IA.2) . As a consequence, this term remains independent of the 
qubit-oscillator detuning A, for which reason its contribution to the {a + a^)^-terms is finite. 

For further convenience, we introduce transformed creation and annihilation operators 
that describe the oscillator-qubit system in the adiabatic limit i2 ^ o, 

.4y|(„+„t,^>y|,„-„t,, (A.0, 

and accordingly, such that [d,a^] = 1. The effective oscillator frequency 



4co 

+ — (A.ll) 

accounts for all quadratic corrections to the oscillator potential in the effective 
Hamiltonian ( |A.9t in terms of the effective operator-valued coupling frequency 

a) = i(AAi-EA|)(j, + ii2An(AA + AE)o:,+§2(cos0c7,-sin0ff^) . (A.12) 

Thus, the effective Hamiltonian can be rewritten as 

J^Q = ri£l{a^d+]^+]^rm^),a,. (A.13) 

Put differently, the qubit-oscillator coupling has been shifted to the effective operator-valued 
oscillator frequency Ci, which depends on the qubit state. 

In order to move fully to the dispersive picture, we also have to transform the system 
operator Q = a + by which the oscillator couples to the environment. Transformation with 
the operator (I A. 2) yields in first dispersive order the position operator 

Q = ''^''Q'^ = {d + d^)-{X^- Ai:)c7, + IX^O- , (A. 14) 

where we have assumed ii « £2. The full system-bath Hamiltonian in the dispersive picture 
finally reads as 

.J^^.M + fiQY,Cn{b„ + bl)+Y,ri(an(bX + ^) . (A.15) 

n n ^ 

Appendix B. Input-output formalism 

In order to compute the response of the oscillator to the external driving, we employ the 
input-output formalism [|20l . which is most conveniently obtained from the quantum Langevin 
equation of the central system Il44lff7l . We derive it from the system-bath Hamiltonian ^ 
via the Heisenberg equation of motion for the bath oscillator coordinates qn= + b^, 

qn + C0^q„^2c„(OnQ . (B.l) 

Here, the oscillator coordinate Q = a + a' enters as inhomogeneity. The formal solution of 
equation (IB. lb for initial time to is 

q„ {t)^qnito) cos{ 0),, (r - fo ) } + sin{ co„ (f - fo) } 

(B.2) 

+ 2c„ / df'sin{co„(f-f')}G(0 , 

JtQ 
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with Pn{tQ) — qn{to)- Inserting this solution into the Heisenberg equation of motion for Q 
yields 

~D?Q-mY^cl f dt'sm{(o„it~t')}Q{t') 

-2nrc„ ( ^„(fo)cos{to„(f-fo)} + ^^^sin{a)„(f-fo)} ) • 
TV ®n / 

For the sake of a compact notation, we define the operator for the incoming fluctuations, 

C(0 =E^" (qn{to)cos{co„{t-to)} + ^^sm{(o„it-to)}) , (B.4) 

which only depends on the environmental operators at initial time and, thus, is independent of 
the central quantum system. 

We replace the sum ^„ |c„p by an integral over the spectral density J{(o)/n, which for 
the ohmic J{(o) = aco becomes the derivative of the delta function 5(f — f'), such that the 
time integral can be evaluated. In doing so, we arrive at the quantum Langevin equation 

Q + 2a^lQ + Q."Q = -2a^^7(f), (B.5) 

where we have discarded an initial slip term and a constant potential renormalisation which 
both are not relevant in the present context and beyond transient behaviour Notice that 
dissipation enters via a friction term, while the incoming fluctuations act as stochastic driving 
force. 

The quantum Langevin equation (IB. 51 ) can also be expressed in terms of the outgoing 
fluctuations by solving the equations of motion ( IB. II ) for q„ with boundary condition at a later 
time time ti > t, i.e., by backward propagation. Then one obtains 

In (f ) = (fi ) COS { 0),, (r - f 1 ) } + '-^^^ sin{ co„ (f - f 1 ) } 

(B.6) 

+ 2c„y^ dt' sm{co„{t-t')}Q{t') . 
The corresponding environment operators define the outgoing fluctuations 

&)=L'^n L(ri)cos{«„(f-r,)} + ^sin{a)„(f-fi)}) . (B.7) 
TV / 

In contrast to (f ), this noise operator depends on the time evolution of the system at earlier 
times f < f 1 . The resulting Langevin equation for the oscillator coordinate Q, 

is characterised by negative damping and the outgoing noise. The difference of both Langevin 
equations links the noise terms via twice the dissipative term by means of the input-output 
relation, which in the stationary limit fo — —°o and f i — > oo reads ll20l 

Ct(0-C(0=2ae. (B.9) 

Even though we have written this relation for a harmonic oscillator, the derivation does not 
rely on particular properties of this system. Thus, equation ( |B.9t is valid as well for non-linear 
quantum systems coupled to an environment. 

If a bath mode is coherently excited by an external driving field, the incoming fluctuations 
are augmented by a deterministic contribution, +Xi^yg{t). Then the input-output 

relation allows one to compute both the averaged outgoing signal as well as its fluctuations 
and noise spectra. 
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Appendix C. Bloch-Redfield master equation 

The numerical data presented in section |4] have been computed with a quantum master 
equation of the Bloch-Redfield type BSl . 

Po(0 = -^k4),Po(f)]-[G,[G,Po(f)]]+ia[e,{e,Po(f)}] , (C.l) 

where 

Q=^J^ dT^ d»cocoth(^^)cos(coT)(2(-T) . (C.2) 

It describes the time-evolution of the reduced density operator po(f) of the qubit plus the 
oscillator The dissipative terms have been derived under the assumption that the bath 
couples weakly to a system operator Q with a vanishing equilibrium expectation value. 
The environment is in a thermal state at temperature T, and the system-bath interaction 
possesses the ohmic spectral density J{(0) = aco with the dimensionless damping strength a. 
Furthermore, X{t) = '^Q{t,tQ)X'W(){t,to) refers to the time evolution of the system operator 
X in an interaction picture described by the propagator '%(f,fo) = exp{iJ^{t — fo)/^}, and 
2 is a shorthand notation for the Heisenberg time derivative i[J%,Q]/h of the system-bath 
coupling operator Q. 
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